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Abstract. We interpret as shear viscosity the anisotropic pressure that emerges 
in inhomogeneous spherically symmetric spacetimes described by the Lemaltre- 
Tolman-Bondi (LTB) metric in a comoving frame. By assuming that local 
isotropic pressure and energy density satisfy a generic ideal gas equation of state, 
we reduce the field equations to a set of evolution equations based on auxiliary 
quasi-local variables. We examine the transport equation of shear viscosity 
from Extended Irreversible Thermodynamics and use a numerical solution of the 
evolution equations to obtain the relaxation times for the full and "truncated" 
versions. Considering a gas of cold dark matter WIMPS after its decoupling 
from the cosmic fluid, we show that the relaxation times for the general equation 
are qualitatively analogous to collision times, while the truncated version is 
inadequate to describe transient phenomena of transition to equilibrium. 



PACS numbers: 98.80.-k, 04.20.-q, 95.36. +x, 95.35. +d 

1. Introduction. 

It is a well known fact that dissipative effects in the context of General Relativity 
must comply with causality and stability requirements [U [51 |31 H]. Also, there is 
an evident theoretical connection between dissipative phenomena and anisotropy or 
inhomogeneity of self-gravitating sources. This emerges from the fact that heat flux 
and shear viscosity couple with the 4-acceleration, shear and spacelike gradients in 
their corresponding evolution (or transport) equations. Since bulk viscosity is the 
only dissipative stress compatible with global isotropy and homogeneity, most articles 
on dissipative cosmological sources deal with the effects of this stress in a Friedman- 
Lemaitre-Robertson- Walker (FLRW) context [5]. However, the literature contains 
also a large number of studies of dissipative cosmological sources under anisotropic and 
inhomogeneous conditions, using Bianchi or Kantowski-Sachs models involving 
heat flux |7j or shear viscosity with the Lemaitre-Tolman-Bondi metric [HI [SI [ID] (see 
also [TT] for inhomogeneous spacetimes with dissipative sources). 

Besides mathematical simplicity, the main justiflcation for preferring a FLRW 
framework, or linear perturbations on a FLRW background, in cosmological studies 
is the conjecture (supported by observations) that the universe is approximately 
FLRW at a large "homogeneity" scale of 150-300 Mpc [12]. Thermal dissipation 
might play a minor role in these scales, as observations seem to reveal that cosmic 
dynamics is presently dominated by (apparently) non-thermal sources (cold dark 
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matter and dark energy [12^). While dissipative phenomena of a thermal nature 
are relevant for understanding early universe interactions, inhomogeneity can be 
safely assumed to be very small in these conditions. Dissipative phenomena also 
arise in self-gravitating (and inhomogeneous) sources at local scales, either stellar 
or galactic [T]. In some cases (inter-stelar or inter-galactic clouds of ionized gas), 
characteristic velocities and energies are basically non-relativistic, but in other cases 
(gas accretion to compact objects or AGN's, jets, photon or neutrino transport) we 
can have non-trivial relativistic and ultra-relativistic effects in conditions of non-linear 
inhomogeneity [T5]. However, as long as we ignore the fundamental physics of dark 
matter and dark energy, we can still try to probe theoretically the possibility of some 
forms of thermal dissipation in these sources and/or their interactions at the cosmic 
scale. 

Fully general inhomogeneity requires numerical codes of high complexity, hence 
we offer in this article a compromise by looking at dissipative phenomena in spherically 
symmetric sources, which in spite of their obviously idealized nature, are still useful to 
examine non-linear phenomena that cannot be studied in a FLRW framework or with 
linear perturbations. By considering "LTB spacetimes" that generalize to nonzero 
pressure the well known LTB dust solutions [2], we obtain a class of spacetimes 
that can be fully described by autonomous first order evolution equations that can 
be well handled by simple numerical methods. These models are quite general and 
readily allow for an inhomogeneous generalization of a large number of known FLRW 
solutions. The reader can consult |14j for an extended and comprehensive discussion 
on these spacetimes and their physical and geometric properties. 

The plan of the article is as follows. We describe in section 2 the basic features 
of LTB spacetimes [TTJ [TJ , in which the anisotropic pressure is considered as a shear 
viscous stress [S1IH[TU]. Assuming a conserved particle current and the entropy current 
associated with Extended Irreversible Thermodynamics, we derive in section 3 the 
full causal transport equation for shear viscosity [TJ [21 [3l H] . In section 4, we provide 
the "fluid flow" evolution equations for LTB spacetimes [TS], equivalent to the field 
equations, in terms of suitably defined quasi-local variables [T?]. In this description, 
the local thermodynamical state variables are gauge invariant "exact" perturbations 
of their quasi-local equivalents. In section 5 we specialize the evolution equations for 
a local equation of state corresponding to a generic ideal gas that covers the cases 
of (i) a classical ideal gas and (ii) the coupled mixture of a non-relativistic gas and 
radiation (the "radiative gas" [TJ [21 [HI HH] ) ■ In section 6 we specialize the evolution and 
transport equations for the ideal g model of a gas of non-relativistic WIMPS 

after their decoupling from the cosmic fluid [H [12 US] , when particle numbers are 
conserved. We evaluate the relaxation times for the full transport equation and for 
its "truncated" version (the Maxwell-Cattaneo equation). In section 7 we compare 
numerically these times with mean collision times, showing that they are qualitatively 
analogous in the relaxation time scale. These numerical examples also show that the 
truncated equation is inadequate to describe the transient phenomena of transition 
to equilibrium for gas of WIMPS. This result is analogous to that obtained for the 
decoupling of matter and radiation in the radiative gas [10] . We summarize the results 
obtained in section 8. 



Shear viscosity, relaxation and collision times in spherically symmetric spacetimes. 3 
2. LTB spacetimes in the "fluid flow" description. 

Spherically symmetric inhomogeneous dust sources are usually described by the well 
known Lemaitre-Tolman-Bondi metric [TTJ [T3] 



ds^ 



-dr 



f + sin^ 



(1) 



1-K 

where R = R{ct,r), R' — dR/dr and K = K{r). A large class of spherically 
symmetric spacetimes follow at once by considering the most general source for m 
a comoving frame (u'^ — Sq), which is the energy-momentum tensor 

T'^'' = + + (2) 

where fi and p are the matter-energy density and the isotropic pressure, h""^ = 
the induced metric of hypersurfaces orthogonal to u° , and 11°'' is the 
symmetric traceless tensor of anisotropic pressure. We will call "LTB' spacetimes" to 
all solutions of Einstein's equations for ([T|) and ([2| . 

Besides the scalars /i and p, and the tensor 11°^, the remaining basic covariant 
objects of LTB spacetimes are: 



e = v„u" = 



2R R' 

ll^R'' 



2{KR)' 



<Jnh = V 
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ah 



UcUdC 



abed 



Expansion scalar, 

Ricci scalar of the ^T, 

Shear tensor. 
Electric Weyl tensor, 



(3) 

(4) 

(5) 
(6) 



where R = u°-\7aR, Vq = /i^Vb, and C'""^ is the Weyl tensor. 

For spherically symmetric spacetimes, the symmetric traceless tensors ct"'', H"^ 
and E""^ can be expressed in terms of single scalar functions as 

where = — 'irj'^rf' and 77" = VhT^S'^ is the unit vector orthogonal to m° and to 
the 2-spheres orbits of S0(3) parametrized by (0,0). 

The field equations G"'' = kT"'' (with k = SttG/c^) for Q and ^ are 

R{if + K)] 
1 

" 3 



K n R^R' 
npR^R' 
R' 



R 
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From (|3|, ([8a|-(8cl and ^ we obtain the expressions for E and S in terms of metric 
functions 
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R R' 
R~ R/ 
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-V u 
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R^ + K 
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The energy-momentum balance equations VbT"'' = for ^ are 

(/i + p)e- dabn"' = -(^i + p) e - 6 s p, 



^aP 



p' _ 2V' = 6V 



R^ 



(9) 

(10a) 
(106) 
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so that pressure gradients are effectively supported by the anisotropic pressure. 

Bearing in mind Q and the remaining previous equations, all covariant objects 
(scalars and proper tensors) in LTB spacetimes can be fully characterized by the 
following set of local covariant scalars: 

{/i, p, T', e, E, f , 37^}. (11) 

Given the covariant "1+3" time slicing afforded by the evolution of these scalars 

can be completely determined by the following set of "fluid flow" scalar evolution 
equations |T5] 

e- -^-^(/i + 3p)-6S2, (12a) 

= _ {^i+p)e (12&) 

±= -^S + E^-f+^P, (12c) 

e 



(12d) 



together with the spacelike constraints 

IT 
R 

3 / ■ " " i? 



{p-2V)' ^6V^ = 0, (13a) 

(r)\' fi" 

^ 3S^ = 0, (13&) 



fi+lv] +£' + 3£^ = 0, (13c) 



K 

6 V ' 2' J ' ' ' ^' R 



and the Friedman equation (or "Hamiltonian" constraint) 



3j 3^-ir+^' 

The system ( 12a|-( p^ is e quivalent to the field plus conservation equations VbT"^ — 
(equations ( 126l and ( 13al). However, this system requires an equation of state linking 
/Li, p and V to become determined, and the time and radial derivatives (in general) 
do not decouple. Hence, we will consider in section 5 another set of equivalent (but 
easier to handle) scalar evolution equations. 



3. Extended Irreversible Thermodynamics. 

In order to arrive to a determined set of evolution equations, we need to prescribe an 
equation of state that is suitable for a given physical model. If the desired model is 
a thermal system, it is necessary to consider /i and p as thermodynamical scalars. 
In particular, a very useful system is the ideal gas associated with the following 
equilibrium equation of state [1] [S] 

lj, = mc^n-\ kT = — (15) 

7 — 1 n 

where n is the particle number density for a gas whose particles have mass to, T is 
the temperature, k is Boltzmann's constant and 7 is a constant. 
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/i = rnc n+ -p, 



For 7 — 5/3, the generic equation of state (151 becomes [8l [9l [10] 

^P: kT=- 
2 n 

which is the equation of state of a non-relativistic hmit of the classical ideal 
gas (Maxwell-Boltzmann gas). Another system that can be described by (15) is 



(16) 



a suitable approximation to a mixture a non-relativistic and an ultra-relativistic 

gas [numiaiioj: 



kTi 



TO(nr)'T'(i, 
Pin 



3p(ur), 



(„r) 



kT, 



n 



(nr) 



(ur) 



P(ur) 
'^(ur) 



(17) 

where subindices and (m.) respectively stand for non-relativistic and ultra- 
relativistic. If we assume that P(nr) P(ur)j but the non-relativistic gas provides 
the major contribution to rest mass (m(„r)n(nr) ^ W(ur)?^(ur)), then (17) becomes 

P(ur) 



Hnr) 



which is the equation of state (15) with 7 — 4/3 



"(ur) 

In practice, one uses ([T8| 



(18) 



to 



describe the so-called "radiative gas" , which is a tightly coupled mixture of baryons 
and photons described as a single "dust plus radiation" fluid. In particular, since we 
neglect thermal motions of non-relativistic particles, ?7i(ni.)C^ ?^(nr) could be the rest 
mass density of cold or "warm" dark matter and non-relativistic baryons, and so TO^^r) 
could be taken as the mass of a neutralino or another supersymmetric DM particle 
candidate. 



For either form (16) or (18), we will assume particle number conservation 



(19) 

hence, if we consider a dark or warm DM gas described by ( 16 1, we would be necessarily 



looking at dissipative effects after the "freeze out" era, when thermal equilibrium is no 
longer kept by particle annihilation [T^[T31[T2]. On the other hand, considering the 
radiative gas model, then (18) with particle conservation is appropriate to describe 



the photon-electron interaction associated with Thomson or Compton scattering. 

Since (p| ) contains anisotropic pressure, which is not involved in the equation 
of state (15), it is natural to consider this pressure as a shear viscous stress 



associated to irreversible processes, whether in the classical ideal gas of WIMPS 
(16 1 or in the radiative gas (18). Considering the fact that Extended Irreversible 



Thermodynamics (EIT) provides the most advanced theory complying with causality 
and stability [U El [3 H] , we construct an entropy current 5*° within the framework of 
this theory. Since = for LTB spacetimes and the only dissipative stress is shear 
viscosity, the entropy current is 



^(cq) 



5'(cq) 



r\nT 



(20) 



2r]nT 

where we used ([7|, and r, r\ are, respectively, the relaxation time and the coefficient 
of shear viscosity, while the specific entropy, S'^'^'^ , is given by the equiUbrium Gibbs 
equation 



rd5(<=^) = d(^)+pd('^], 



(21) 
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so its projection with respect to u° and the balance equation (10a I yield 



(22) 



The condition VaS"" > 0, together with (22 1, leads to the transport equation for shear 
viscosity [TJ lU d U 



CThlhfilcd 



2t]T 



2riaab = 0, 



(23) 



where eg = 0, 1 is a "switch", so that (23 1 is the complete transport equation if eo = 1, 
and we get the "truncated" or Maxwell-Cattaneo equation if eq = 0. Using ([3| and 
^ , equation ( 23 ) becomes the following scalar equation 



ctV + V + 2t]J: 



2 



T 

f 



e 



= 0. 



(24) 



To apply EIT to the non-relativistic and radiative gases, we need to substitute 
the equation of state ( [l5| and utilize the forms of the coefficient of shear viscosity for 
these gases. From [T][2n3lllj, we have 

1, non-relativistic ideal gas 
|, radiative gas 



77 = apcT, 



(25) 



Hence, inserting p — nkT and the particle conservation law (19 1, the transport 



equation (24 1 becomes 



r + 2apY, + eoV 



r = o, 



(26) 



which clearly reveals how the difference between the complete and truncated equations 
can be dynamically significant, as it involves the term V T/T. The entropy production 



subjected to the conservation law (19 1 follows readily from (20 1 and (22) as 

752 



VaS"" = nS = 3kn 



CT p 



ap-^ 



(27) 



where we used ( 25 ) and ( 26 1 to eliminate V 



In order to examine (20) and (26) we need to solve the field equations, or 



their equivalent "fluid flow" evolutions equations (12al-(14l, which would render the 



functional forms of the involved thermodynamical scalars. We look at this matter in 
the following section. 



4. Quasi local evolution equations. 



We can obtain an alternative set to the evolution equations (12a)-(14l that is 
completely equivalent, but easier to deal with numerically |14j . This follows from 
using instead of the local scalars (11 1, the scalar representation given by quasi-local 
variables defined by the map 

/o AR'^R'dx 



J, ■.X{V)^X{V), 



(28) 



R^R'dx ■ 

where X{'D) is the set of all smooth integrable scalar functions A defined in any 
spherical comoving region D of the hypersurfaces '^T orthogonal to u", containing a 
symmetry center marked by r = 0. The functions A* : I? — > R that are images of J7* 
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will be denoted by "quasi-local" (QL) scalars. In particular, we will call A^, the QL 
dual of A. See [T3] for a comprehensive discussion of of the map (pSl). 

Applying the map ( 28 1 to the scalars Q and ^TZ in ^ and we obtain their 
QL duals 



3R 



6K 



(29) 



Applying now ( 28 1 to fi and p, comparing with 



field equations transform into 

2 



-{8b \, and using (29 1, these two 



(30a) 



02 ^ 



(306) 



which are identical to the FLRW Friedman and Raychaudhuri equations, but among 
QL scalars. These equations can be further combined to yield identically the FLRW 
energy balance equation: 

/i* = - (fi* +P*) 6*. (31) 

so that the QL scalars {/i*, p*, <^^.} effectively satisfy FLRW evolution laws. 

In order to relate local scalars to their and QL duals, we introduce the following 
"relative deviations" or "perturbations" 

A-A^ 



<5(^) = 



A = A, 



1 + <5(-4) 



(32) 



Therefore, all scalars A in (111 can be expressed in terms of their duals A* and 
perturbations (S'"*); 



M = A** 



1 + ^(p) 



e = 



= ^7^* 



whereas E, V and £ follow as 
1 



V 



[6 - e,] 

o 

2 [P-P*]=2 



3 

P. s^p\ 



£ = - 



6 



^P 



(33) 

(34a) 
(34&) 
(34c) 



which leads to an alternative QL scalar representation {A^, (5'"^'} that it is fully 
equivalent to the local representation. We derive now the evolution and constraint 
equations for this representation. 

From differentiating both sides of ( 28 1 and using ( 32 1 , we can relate radial 
gradients of /i* , p^ and with their corresponding S functions by 



a; 



3R' 



3i?' 



P*^ 
p* 



3i?' 



(35) 



while (30&I and (311 are evolution equations for /t* and 0*. Hence, the evolution 
equations for (5^^^ and S'^^^ follow from the consistency condition [A^,]' apphed 
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to (30&I, (31 1 and (351 for A = Q^, and fi^,. The result is the following set of 
autonomous evolution equations for the QL scalar representation {A^, (5''^^}: 



= -[n 
_ el 

3 



^ [1 + 3u;] Ai*, 



(36 a) 
(36&) 
(36 c) 



e. 



6(e./3) 



's(e)_s(p) + 3w(^S<^''^-d^P^)], (36d) 



where w; = 

The constraints associated with these evolution equations are simply the spatial 



gradients (351, while the Friedman equation (or Hamiltonian constraint) is (30a 



Notice that ( 35 ) follow directly from differentiating the integral definition ( 28 1 , so by 



using the QL variables we do not need to solve these constraints in order to integrate 



(36al-(36dl 



It is straightforwards to prove that the evolution equations (36a|-(36rf) and the 



constraints (30a) and (35 1 are wholly equivalent to the evolution equations (12a)-(14l 
of the fluid flow formalism of Ellis, Bruni and Dunsbury [15]. It is also important 
to mention that the QL representation {A^., i5*-^-'} leads to a characterization of LTB 
spacetimes as exact, non-linear, gauge invariant and covariant perturbations over 
a FLRW formal background defined by the QL scalars A^,, which satisfy FLRW 
dynamics. See [T4j for details. 



5. Evolution equations for the generic ideal gas. 



In order to integrate the system (36al~(36dl we need to prescribe a relation between 
p* and S^^\ S^'P\ Since we are interested in thermal dissipative phenomena 
characteristic of a hydrodynamical regime of short range interactions, the physically 



meaning full equation of state ( 15 ) is that relating local variables n and p, and not 
QL variables. However, (151 is a linear functional relation, hence its validity as a 



local relation and the assumption of particle current conservation ( 19 1 are sufficient 



conditions to render (36a)-(36dl a fully determined system in which the QL variables 



are basically auxiliary variables (and the physical variables are the local ones). 

Assuming the local equation of state (15 1 and using (32 1 with A = ji, n, p leads 



to the following conditions on the QL variables 

P* 



m (? n» 



7- 



1' 



(37a) 



(7-1)a^. 



. j(p) 



(376) 

Using the particle numbers conservation law ( 19 1 with n ^ (l + ^(")), together with 



(37al-(37&), transforms (|36a 

= - 7 p* 6 

3 2 ^' 



{36 d\ into the fully determined system 



= - 



71 P* 



(38a) 
(38&) 

(38 c) 

(38 d) 
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sip) 



mc^ n 



7i 



(38e) 



(38i) 



with 7i = (37 — 2)/(7 — 1). Once the system (38aH(38/l is solved numerically for 
appropriate initial conditions (see appendices of [14 ), we obtain the local variables 
p, O, S, V from (33 1, (34a) and (34&I, while the temperature T follows readily as 



1 + 



(39) 



With the help from (33), (34al, (34&), (38a|-(38yi and (39), the transport equation 



( 26 1 reduces to the following two algebraic constraints defining the relaxation times 
1) and truncated (e = 0) cases: 



for the full (e 



cr : 



e*(5(e) [4a((5(p))2 + (3 + 8a)S(P) + 37 + 4a] ' 
3(5(P) 

CT — 

[((4a + 3)(5(e) + 37)5(P) + (4a + 3j)S(^)] 
while the entropy production law (p7| leads to 



3/c((5(P))2 
4a [1 + 



1 

CT 



(ep - 1)94(7 + + (1 + (5(Q))7' 

1 + 



= 1, (40) 
e = 0, (41) 

(42) 



However, substituting r from either (40) or (41) into (42) we obtain the same 

37] 



expression of S for the full and truncated cases: 

fee* (5(®) 6(-P^ [4a {S^P^f + (3 + 8a) S^p^ 



S- 



4a 



(43) 



4a [1 + (5(P)]3 

Dissipative effects associated with shear viscosity for thermal systems associated with 



(15) can be now examined by usin g th e numerical solution of (38aH(38J) to calculate 



from (42) 



the relaxation time scale given by (40) or (41), as well as the entropy production nS 



6. The gas of WIMPS 

Dissipative phenomena associated with shear viscosity in spacetimes with LTB metrics 
have been studied mostly on the radiative gas model [S]|TU] (but see [5]). In particular, 
the comparison between relaxation and collision times was examined in [T^ for this 
model in the context of the cosmological radiative era. In this article we consider 
the same issue, but for a gas of non-relativistic cold dark matter particles (WIMPS) 
after its decoupling or "freeze out" from the cosmic fluid, when thermal equilibrium 
is no longer maintained by particle annihilation [T5J [T31 [TS] . Since cold dark matter 
has no effect on cosmic nucleosynthesis, this decoupling must have happened before 
nucleosynthesis at around t ^ 200 sec, and so the gas of WIMPS can be described 
as an ideal gas in the non-relativistic limit, corresponding to the equation of state 
( 16 1. Hence, the expressions for the coefficient of shear viscosity, relaxation times and 
entropy production are (25), (37a)-(43l, for the the values 7 = 5/3 and a 



1. 



The relaxation time is a mesoscopic quantity that could be, in principle, obtained 
by means of collision integrals but cannot be given in terms of an "equation of 
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Figure 1. Coldn ess parameter 13. The figure displays the function /3 = 
ra(? l(kT) given by ( |5l[ l for the ideal gas of WIMPS configuration described in 
section It] The layers near the center (r = 0) bounce and collapse to a black hole 
where thermal motions dominate rest mass (/3 — > 0), though the hydrodynamical 
regime is no longer valid in this stage. 




Figure 2. Relaxation vs Hubble times. The figure depicts the logarithm of 
T for the full transport equation (A), for the truncated version (B), the Hubble 
time 3/© (C) and 5 (D). The vertical dotted line depicts the extension of the 
relaxation time scale up to t ~ 3/©. Notice how for central (over-density) layers 
(left panel), with more thermal energy (lesser /3), this time scale has a much larger 
extent than in the layers corresponding to the cosmic background (right panel). 
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Figure 3. Relaxation vs collision times. The figure depicts the logarithm of 
T for the full transport equation (A), the Hubble time 3/0 (C) and collision times 
(dotted curves) given by l |46| with the numbers indicating the exponent sq, and 
with the WIMP mass— energy of 100 GeV. Notice how t is qualitatively analogous 
to collision times with cross section areas a ~ 10~^®cm'^ consistent with weak 
interactions of cold dark matter WIMPS. As in figure 2, the relaxation time scale 
has a larger extent in layers in the over-density (left panel) than in the cosmic 
background (right panel). 



state" relating macroscopic thermodynamical scalars. Usually, this quantity is taken 
simply as a mean collision time, or it is assumed to have the same order of magnitude 
value as these times. However, as shown by the results of [TU] in the radiative gas 
model, there is no reason for this to be the case. Since these two time scales follow 
from physically distinct concepts, they must be different functions that could exhibit 
qualitatively analogous behavior and/or could be of the same order of magnitude. 
For an ideal gas the mean collision time is given as [TJ [131 dS] 

= n\xM\ ^ (44) 

an an* (1 + d*."^) 

where a = a{n, T) is the collision cross section area, whose precise functional form 
follows from the specific particle interactions involved in the gas model. For a gas of 
WIMPS, we can identify a decoupling stage as cosmic times ct = ct-o for which the 
reaction times compare with the Hubble expansion time tn 

3 

na{n,T) ^ ctH ^ —, (45) 

so that for t < tiy, before its decoupling from the cosmic fluid, a is associated 
with particle pair annihilations and its form follows from theoretical considerations 
pertinent to supersymmetric cold dark matter candidate particles |13l Moreover, 
we will examine dissipative effects in the gas of WIMPS for t > t^, after this freeze 
out when particle numbers are conserved. The justification for these after freeze 
out dissipative processes comes from the assumption that there could have been 
dissipation in the earlier stage t < tj^, and so it is reasonable to assume that once 
particle annihilations stop at ( 45 1 , there should be a short timed relaxation process 
characterized by a weak self-interaction associated with very small cross section 
areas, so that after this process the fluid becomes completely non-coUisional. Since 
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this relaxation should be of short duration, we can model these cross section areas 
empirically by the simple ansatz |16l 113] 

-40 < So < -34. 



10"° cm^, 



(46) 



Hence, 



we expect r in (40) and (41 1 to exhibit a qualitatively analogous behavior 
as (44) for cross section areas having magnitudes given by (46). In particular, the 



existence of an interaction that can be associated with shear viscosity requires that 
these time scales are of lesser magnitude than the Hubble time: ct < 3/0 and 
cicoi < 3/0, with the relaxation time scale given by the cosmic time ct such that 
cr ~ 3/0 and ctcoi ~ 3/0, and thus, cr ^ cicoi at these cosmic time values. Hence, 
when CT > 3/0 the gas expands in a non-coUisional stage. However, for earlier times 
CT and ctcoi need not be the same function, just have comparable magnitudes. Also, 
for the relaxation time scale in which ct < 3/0, we must have necessarily S > 0, so 
that there is entropy production with —> as cr and c^coi overtake 3/0 and entropy 
becomes a maximum associated with equilibrium conditions. 

In order to test numerically these conditions, we define the following dimensionless 



variables associated with n^, and O* in (38al~(38yi (notice that the S functions are 
already dimensionless): 



Km c 



where Hi ~ l/{cti) is taken as the Hubble scale factor for the initial time surface 
ti ~ 200 sec, and we will consider mc^ = 100 GeV to be the rest mass-energy of the 
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(47) 



WIMP. In terms of (47 1, the collision time (44 1 is given by 

2 9 9 

mc cm 



= 4.3 X 10" 



1 



CTx(l + (5(")) GeY lO'o x{l + 6('''>)' ^'^^^ 

We will examine in the following section these different time scales associated with 



the relaxation scale using the numeric solutions of (38aH(38/l 



7. Relaxation time scales: numeric results. 



In order to set up appropriate initial conditions for x, y and z, we use equations (30a 



and (37a) for 7 = 5/3, leading to 

z1{r) ^ Xi{r) + yi{r) - ki{r), ki{r) 



(49) 



where the subindex j denotes evaluation ai t — ti . Initial conditions for a central over- 
density with small positive curvature that smoothy blends to a cosmic background 
with small negative spatial curvature can be achieved by choosing ki (r) as any smooth 
function for which fci(O) = 0.1 and ki{r) — *■ —0.1 for r —^ oo. Central and asymptotic 
values for Xi and yi are given by 

x,(0) = 1.5, x,(oo) = 0.9, y,(0) = 0.08, 2/,(oo) = 0.02. (50) 



The form of Zi follows from ( 49 1 and ( 50 1 , while the forms for the initi al va lue functions 
[i5'^"']i, [S'^P^i and can be obtained from Xi, yi, Zi by means of (35) evaluated at 

t — ti (see the appendices of [H]). 

An important parameter in thermal systems is the "coldness" parameter 



mc^ _ x[l + S^"'>] 
~kT ~ y [1 + <5(P)] ' 



(51) 
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which, with the numeric vahies of (50l, takes initial values /3i(0) ~ 20 and (3i{oo) ~ 40, 
which are reasonable values for cold dark matter WIMPS that are non-relativistic 
when they decouple at t = to [HI HI US] ■ We display in figure 1 the function (3 
that results from the numeric solution of (38a|-(38J) for the configuration outlined 
above. As the configuration expands we can see how (3 increases for all r to clear 
non-relativistic values 3> 1, but layers in the over-density region (around r = 0) 
collapse to a black hole at around Hict ~ 150, with /3 ^ 0, indicating dominance 
of internal energy density over rest mass energy density near the final collapse. 
However, the hydrodynamical regime is no longer a valid approximation in this stage, 
as WIMP configurations do not evolve to black holes. In more realistic structure 
formation scenarios the WIMP gas becomes non-coUisional and undergoes non- 
coUisional relaxation phenomena, such as virialization [13', leading to stable bound 
structures. 

The relaxation of a viscous dissipative stress requires that S" > while cr < 3/6, 
but both r and tco\ must overtake 3/0 as ^ 0. We test these conditions numerically 
in figure 2, for two different values of r (at the over-density in the left panel and at 
the cosmic background in the right panel) , and for the relaxation times of the full ( 40 1 
and truncated (41 1 transport equations and for S given by (43). As shown by this 
figure, we have ct < 3/0 for all times for the relaxation time (4lJ of the truncated 
equation. Therefore, the relaxation time for the truncated (Maxwell-Cattaneo) does 
not exhibit the appropriate behavior of a relaxation parameter, which means that the 
full transport equation is needed to provide an adequate description of the transient 
dissipative phenomena for the ideal gas of WIMPS. The same result was obtained for 
the radiative gas model in [TU] . On the other hand, the relaxation time ( 40 1 of the 
full transport equation exhibits the expected behavior and overtakes the Hubble time 
3/0 as S ^ 0. We show in figure 3 how the relaxation time (40 1 of the full transport 
equation (for mc^ = 100 GeV) in the whole relaxation time scale is qualitatively 
analogous to collision times with cross sections given by ( [46| with sq ~ —36, which 
characterize expected weak interactions for decoupled WIMPS [TH[T31[TB]. 



8. Conclusion 



We have examined causal dissipation from shear viscosity in the context of a large class 
of inhomogeneous spherically symmetric spacetimes described by the LTB metric ([T]) . 
A generic equation of state was suggested, which contains as particular cases the 
classical, non-relativistic, ideal gas, as well as the radiative gas in the approximation 
in which thermal motions of the non-relativistic species are ignored. We obtained 
a set of evolution equations equivalent to the field and balance equations, whose 
numeric solutions can be used to compute the relaxation times for the full and 
truncated transport equations, the rate of change of specific entropy and collision 
times for suitable cross section areas. We considered the non-relativistic ideal gas as 
an appropriate equation of state for a gas of cold dark matter WIMPS undergoing a 
transition to equilibrium soon after their freeze out and decoupling from the cosmic 
fluid at the outset of cosmic nucleosynthesis. The comparison between relaxation and 
collision times yielded similar results as those obtained in [TD] with the radiative gas 
model, namely, that only the relaxation time from the full transport equation exhibits 
the expected behavior of a relaxation parameter, being also qualitatively analogous 
and of the same order of magnitude as collision times with reasonable cross sections 
for a gas of WIMPS. This result is shown in figures 2 and 3. 
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It is evident that the study of shear viscosity without other dissipative fluxes (heat 
flux and bulk viscosity) is an ideaHzed situation which follows from the constraints 
of the LTB metric. Although the inhomogeneous conditions provided by this metric 
are mathematically tractable, they are not trivial and contain enough structure to 
examine non-linear effects that cannot be studied in a FLRW context or with linear 
perturbation. Another shortcoming is the use the transport equation itself to define 
the relaxation times, as it was done in [HlQllTO], instead of using it as a free parameter 
to be specified. The resulting expressions (|40| and (41 1 are, evidently, approximations 



to the actual relaxation times, but this approximation will be reasonable if the obtained 
quantities behave as a relaxation parameters. As shown in section 7 and in figures 2 
and 3, the relaxation time for the full equation does exhibit the expected behavior, and 
so this approximation is reasonable. Future work along these lines would necessarily 
require a more general metric framework and more elaborated numerical methods. 
This work is presently under consideration. 
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